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ABSTRACT 

We measure the evolution of the specific star formation rate (sSFR = SFR/M sta r) between redshift 
4 and 6 to investigate the previous reports of "constant" sSFR at z > 2, which has been in significant 
tension with theoretical expectations. We obtain photometry on a large sample of galaxies at z ~ 4-6 
located in the HUDF and the Early Release Science fields that have high quality optical and IR 
imaging from HST and Spitzer. We have derived stellar masses and star formation rates (SFRs) 
through stellar population modeling of the rest-frame UV and optical spectral energy distributions 
(SEDs). We estimate the dust extinction from the observed UV colors. In the SED fitting process 
we have studied the effects of assuming a star formation history (SFH) both with constant SFR and 
one where the SFR rises exponentially with time. The latter SFH is chosen to match the observed 
evolution of the UV luminosity function. We find that neither the mean SFRs nor the mean stellar 
masses change significantly when the rising SFR (RSF) model is assumed instead of the constant 
SFR model. We estimate the sSFR at z > 4 for galaxies in two stellar mass bins centered at 1 
and 5 x 10 9 Mq. We find that the galaxies in the lower mass bin have very similar sSFRs to the 
more massive ones (within ~ 0.1 dcx). When focusing on galaxies with Af sta r ~ 5 x 10 9 M Q , we find 
that the sSFR evolves weakly with redshift (sSFR(z) oc (1 + z )°- 6±01 Gyr -1 ), consistent with previous 
results and with recent estimates of the sSFR at z ~ 2 — 3 using similar assumptions. We have also 
investigated the impact of optical emission lines on our results. We estimate that the contribution of 
emission lines to the rest-frame optical fluxes is only modest at z ~ 4 and 5 but it could reach ~ 50% 
at z ~ 6. When emission lines of this strength are taken into account, the sSFR shows somewhat 
higher values at high redshifts, according to the relation sSFR(z) oc (1 + z )°- 9±01 Gyr -1 , i.e., ~ 2.1 x 
higher at z ~ 6 than at z ~ 2. However, the observed evolution is substantially weaker than that found 
at z < 2 or that expected from current models (which corresponds to sSFR(z) oc (1 + z) 2 4 Gyr -1 ). 

Subject headings: galaxies: evolution — galaxies: high-rcdshift 



1. INTRODUCTION 



Large samples of Lyman break galaxies (LBGs) have 
allowed for the study of the propertie s of high red- 
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pies are the result of large investments on high quality 
HST data over deep fields like GOODS (|Giavalisco et al.l 
120041) that have rich complementary multi-wavelength 
coverage. Several studies have explored the observed 
UV and optical colors of these galaxies. Through 
the standard technique of synthetic stellar population 
modeling of the observed spectral energy distributions 
(SEDs), the physical properties of these galaxies, such 
as the stellar mass (M star ), star formation rate (SFR), 
dust a ttenuation age, etc, have a lso be en explored 
(e.g.. lYan et al.l 120061: lEyles et all 120051: IStark et all 
2001 IGonzalez et all 120101 120111: IPapovich et all 120111: 



Lee et al .1 120121 : iBouwens et alJl2009l 120121 ). 
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These studies have shown that, at rest-frame UV and 
optical wavelengths where LBGs are more amenable to 
observations, the SEDs show very similar colors with 
weak trends of bluer colors as a functio n of decreasing, 
UV luminosity and increasing r edshift (IBouwens et all 
[20091 120121 I Wilkins et alll20ll iFinkelstein et all 120111 : 
but see also lDunlop et al.ll2012l ). As a direct consequence 
of this, the physical properties estimated through SED 
fitting, using simple models, are also remarkably similar. 

Particularly intriguing are the results that indicate 
that the specific star formation rate (sSFR) of sources 
with a given M star re mains approximate l y constant from 



2 to 



7 (IStark et all 120091 



Gonzalez et al 



120101 : IMcLure et al.ll2011l) . In particular. IGonzalez et al 
(|2010h shows that for sources with M s t a r = 5 x 1O 9 M0, 
the sSFR shows no evidence for significant evolution 
(sSFR ~ 2 Gyr" 1 ) from z ~ 7 to z ~ 2. Such a result 
is at odds with the fairly generic theoretical expectation 
that the sSFR should decrease monotonically with cos- 
mic t ime (e.g., IWeinmann et al.l 1201 lh iKhochfar k. Silkl 

IllD). 

Although suggestive, much of the early work pointing 
toward a relatively constant sSFR was limited in many 
important aspects and based on a number of simplifying 
assumptions. For example, early work assumed zero dust 
extinction for the SFRs of z > 4 galaxies. This was moti- 
vated by the very blue UV colors observed in the galaxy 
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SEDs, which have subsequently been measured more ac- 
curately thanks to the better wavelengt h coverage pro- 
vided by HST /WFC3 observations (e.g., iBouwens et al.l 
I2012T ). Early works also adopted exponentially declining 
or constant star formation histories (SFHs). This is in 
apparent contradiction with the evolution of the UV lu- 
minosity function (LP), which suggests that rising SFHs 
are a better match for the evolution of g alaxies in early 
cosmic times (e.g., iPapovich et al J 120 111 ). Finally, much 
of the early work did not correct for the effect of optical 
emission lines on the estimate of stellar masses (M sta r), 
and hence on the sSFRs. 

Several recent studie s have attempted to r edress some 
of these shortcomings. IBouwens et al.l (|2012D . for exam- 
ple, used new measures of the dust extinction (based 
on UV col ors) to correct the ear l ier sS FR estimates at 
z ~ 4-7. iSchaerer fc de Barrosl (|2010T l considered the 
possible effects of the optical emission lines on the stel- 
lar masses and sSFRs derived for high redshift galax- 
ies, and this effort ha s been considerably extended i n 
Ide Barros et~aTl (|2012l see also iCurtis-Lake et al.ll20l"2t) . 

In this paper we attempt to bring more consistency in 
the exploration of the aforementioned issues, exploring 
the effect of dust reddening in the UV colors, SFRs, and 
M st ar measurements, as well as the effects of choosing a 
SFH that better matches the evolution of the UV LF, 
simultaneously. The goal is to make estimates for the 
physical properties of z > 4 LBG that use empirically 
motivated assumptions that better match a larger range 
of observations. We plan to investigate whether this new 
set of assumptions still shows an approximate plateau 
in the sSFR at z > 2 or shows evidence for an increase 
toward high redshift as predicted by theory. 

We now provide a brief plan for this paper. In Section 
2, we briefly describe the observational data and selec- 
tion criteria used here. In Section 3, we describe our 
approach to stellar population modeling, detailing the 
specific assumptions that we make and the effects these 
assumptioins have on the physical properties we derive 
for LBGs. In Section 4 we present the new measurements 
of the sSFR at high- redshift. We discuss the results in 5, 
and summarize our findings in Section 6. 

Throughout, we use a (Ho, Am, ^a) = (70km s _1 , 
0.3, 0.7) cosmology when n ecessary and we quo te all mag- 
nitudes in the AB system (|Oke fc Gunnlll983l) . 

2. DATA 

This work is based on a large sample of star form- 
ing galaxies in the z ~ 4 — 7 re dshift range found us- 
ing the Lyman Break technique (jSteidel et al.lll996f) in 
the u ltra-deep HUDF/WFC3 field (e.s.. besch et alJ 
20101) and the deep Early Release Science (ERS; 
Windhorst et al.1 1201 If) fields. B oth fields have deep 
HST /ACS (|Giavalisco et"all[200l and WFC3/IR imag- 
ing, as well as the 23 h Spit zer /IRAC d ata from the 
GOODS program (jGiavalisco et al.l l200l. T h e sam - 
ple was already pre sente d in iGonzalez et al.l (|2011h . 
IBouwens et~aTl |2007l ). and IBouwens et al.l (|2012l) but we 
provide a short summary here. 

2.1. HST ACS and WFC3/IR photometry and sample 

selection 

Over the ERS field, both the ACS optical 
(-8435^606*775^850) and the WFC3/IR (V 098 J 125 H 16 o) 



TABLE 1 
Sample Summary. 





2 ~ 4 


z ~ 5 


z ~ 6 


ERS 


270 / 524 


77 / 123 


20 / 32 


HUDF 


137 / 205 


41 / 55 


38 / 51 


TOTAL 


407 / 729 


118 / 178 


58 / 83 



Note. — Number of sources in our z ~ 4, 
z ~ 5, and z ~ 6 samples. IRAC photome- 
try of these sources requires fitting and sub- 
traction of the flux from surrounding fore- 
ground neighbors. This is not possible in 
all the cases. The table shows the number 
of sources in our samples with clean IRAC 
photometry over the total number of sources 
in these samples. 

images from HST reach depths of ~ 28 mag (5a 
measured on 0'.'35 diameter apertures). In the HUDF 
field, the ACS optical data (BViz) are 1.5 — 2 mags 
deeper, and the WFC/IR data (Y105 Ji25-ffi6o) ar e 1.5 
times deeper than the ERS. 

The HST photometry fo r these sources was perf ormed 
using the SExtractor code (jBertin fc Arnoutill996f ). The 
images were registered to a common frame and then PSF- 
matched to the -ffi6o-band data. SExtractor was run in 
double- image mode with the detection image constructed 
from all images with coverage redward of the relevant Ly- 
man break. Fluxes are measured using Kron-style pho- 
tometry. The criteria used to select sources are as fol- 
lows (for details and contamination rates please refer to 
IBouwens et~aTll2007D : 

z ~ 4 B-dropouts: 

(£435 - V 6 06 > LI) A [B435 - V 60 6 > (Veoe - 2 85 o) + LI] 
A (V 60 6 - 2850 < 1.6) 
z ~ 5 V-dropouts: 

{[V 60 6 -«775 > 0.9(^775 -Z 8 50) + L9] V ( V 606 - «775 > 2)} 

A (Veoe - «775 > L2) A {i 775 - z 850 < 1.3) 
z ~ 6 i-dropouts: 

(«775 - ^850 > L3) A (2 8 50 ~ ^125 < 0.8). 

The combined samples contain a total of 729 sources 
at z ~ 4, 178 at z ~ 5, and 83 at z ~ 6 (see table [T]) 

2.2. Spitzer /IRAC photometry 

While the HST /ACS and WFC3/IR data allow us 
to find LBGs at z > 4 and to study their rest-frame 
UV light, data from Spitzer /IRAC is needed to sam- 
ple the rest-frame optical light from such high redshift 
galaxies. These constraints are cr ucial to derive reliable 
stellar masses from SED fitting (jPapovich et al.1 l200ll; 
lLabbe et al1l2010al) . 

The HUDF and the ERS fields were imaged with 
Spitzer in the four IRAC channels as part of the Spitzer 
GOODS program (Dic kinson et al.1 12003:. In this work 
we only make use of the two most sensitive channels cen- 
tered at 3.6 and 4.5 ^tm respectively. The images have 
integrations of ~ 23.3 h (the HUDF was imaged twice, 
thus amounting to a total exposure of ~ 46.6 h). The 
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depths of the images were measured dropping 2'.'5 diam- 
eter apertures in empty regions of the image and deter- 
mining the RMS. After applying a 1.8 x flux correction 
to account for the light outside such aperture in the case 
of a point source, the depths correspond to 27.8 mag (la) 
for a single 23.3 h image in the [3.6] channel and 27.2 mag 
for the [4.5] channel. 

The size of the IRAC PSF, however, is too broad to use 
standard photometric techniques on our sample of faint 
LBGs. The flux within an aperture is expected to be 
contaminated by flux from neighboring sources spilling 
over onto the same aperture. The neighboring sources 
are in general brighter and the contaminating fluxes can 
be significant. The availability of higher resolution im- 
ages (though at different wavelengths) with HST allows 
us to model the source and its neighbors and subtract off 
the expected contaminating flux in our apertures. This 
metho d has already been descri be d in several previous 
work s (lLabbe et alj|2006t l2010allbl; iGonzalez et~aTl l2010l 
[201lHWuvts et al.ll2007Hde Santis et al.ll2007t) . After the 
area around the source has been cleaned from the flux 
of neighbors, we perform standard aperture photometry 
using 2'.'5 diameter apertures. An aperture correction is 
derived from the higher resolution "template image, in 
this case the WFC3/IR -ffi6o image. This image is con- 
volved to the resolution of IRAC and an aperture correc- 
tion factor is estimated for a 2'.'5 diameter aperture. This 
correction factor is then applied to the fluxes measured 
with IRAC and involves multiplying the flux of sources 
(in 2'.'5 diameter apertures) by 1.8 in both the [3.6] and 
[4.5] IRAC channels, consistent with the correction one 
would apply to point sources. 

The cleaning procedure cannot always produce accept- 
able models of the source and neighbors, resulting in poor 
subtractions for approximately 40% of the cases. The 
individual cleaned stamps were inspected by hand and 
those with poor subtractions were excluded from the sub- 
sequent analysis. This criterion mostly depends on the 
feasibility of modeling the neighboring sources, and as a 
consequence, is not expected to introduce any significant 
biases in our results. 

3. STELLAR POPULATION MODELING 

It has been customary to study the physical properties 
of high-rcdshift galaxies by comparing their observed UV 
and optical SEDs with synthetic stellar population (SSP) 
models. The quality of the flux constraints is generally 
not adequate to constrain all model parameters simulta- 
neously. In particular, models with a variety of mctallic- 
ities, initial mass functions, and star formation histories 
can all reproduce the flux constraints almost equally well. 
As a result, these model parameters must be constrained 
independently or fixed to reasonab le values . In th e fol- 
lowing analysis we have assumed a iSalpeterl £l955) IMF 
with cutoffs at 0.1 and 100 Mq. 

It is expected that the metallicities of high-redshift 
galaxies be some what lower than that of local galax- 
ies. For example, iMaiolino et a l. (2008), from a small 
sample of 2 ~ 3.5 galaxies determines gas metallicities 
of ~ 0.2 Zq. A number of other observa tional studies 
have also found similar metallicities (e.g., iPettini et al.l 
l2000l) . and 0.2 Zq is also expected in many carefully- 
constr ucted cosmological hydrodynamical simulations 
(e.g., iFinlator et al1l20lTI ). 



In our analysis we therefore assume a fixed metallicity 
of 0.2Z Q . The resulting SFRs and stellar masses, which 
are the focus of this work, are only weakly sensitive to 
this assumption. 

For the SFHs, most early studies have generally as- 
sumed models with SFRs that decline exponentially as a 
function of time with a characteristic timescale r left 
as a free parameter. The particular case r — > 00, 
which corresponds to models with constant star forma- 
tion (CSF), is sometimes used as the fiducial model. Be- 
cause in a smooth SFH the UV luminosity can be directly 
linked to the SFR, these SFHs would predict a UV LF 
that shifts to brighter luminosities at increasing redshift 
(or one that remains constant in the CSF case), which 
is at odds wi th the current determinations of the UV 
LF at z > 4 (iBouwens et alJl2007t iMcLure et "all 120091: 
lOesch et al.ll2010l : lBouwens et al.ll2012t ). If the SFHs. are 
in fact smooth, then the observed overall dimming of the 
UV LF predicts a me an SFH that is better character- 
ized by a rising SFR (iPapovich et~"aTl 120111: iStark et al.l 
I2009t iReddv et al.ll2012t) . Oesch et al. (in preparation), 
shows that the evolution of the UV LF can be well re- 
produced by an exponentially rising SF H with a charac- 
teristic timesca l e r ~ 500 Myr (see also iSmit et aT1l2012l : 
iPapovich et al.l [20111 ). In deriving the physical proper- 
ties of galaxies in our samples, we will look in detail at 
the extent to which these parameters depend on whether 
we assume a constant or rising star formation (RSF) his- 
tory. We particularly examine the effects on the sSFR 
and its evolution with redshift. 

3.1. Reddening 

In addition to the necessary simplifying assumptions 
described at the beginning of this section, other compli- 
cations remain in the SED fitting process. In particular, 
it is known that the effects of reddening produced by 
dust and the reddening produced as a consequence of 
the aging of the population are largely degenerate. This 
does not have a substantial impact on the stellar masses 
which generally remain constrained within ~ 0.3 dex for 
a range of ages, but the varying dust corrections have a 
direct impact on the determination of the SFRs. 

This degeneracy cannot be broken with the current 
data. However, one potentially fruitful way forward for 
us is to not consider the full range of stellar population 
models available to us (with ranging ages and dust ex- 
tinctions) and to estimate the dust extinction for indi- 
vidual galaxies based on the UV continuu m slope we ob- 
serve using known IRX-/3 relations (e.g., iMeurer et al.l 
Il999f ). Not only does this approach have some justifi- 
cation based on observations of z ~ 0-2 galaxies (e.g., 
IMeurer et al.lll999tlRe"ddv et al.ll2012| ). but following this 
approach we are able to naturally reproduce the UV 
color-luminosity and UV-to-optical color-luminosity re- 
lations. 

Figurcs[T]and[2]show the observed colors of our z ~ 4-6 
samples. The top panel of Figure Q] shows the observed 
!775 — [3-6] color as a function of observed [3.6] mag. The 
different symbols are used to divide the IRAC detected 
(2a 1 solid blue circles), the marginal detections (l-2<r, 
open squares), and the non-detections (open gray cir- 
cles). The histogram shows the color distribution of the 
sample. The two solid horizontal lines show the minimum 
color (corresponding to a > 10 Myr old population) and 
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Fig. 1. — The observed colors of z ~ 4 LBGs. Only sources with 
reliable IRAC photometry are included (sec text, section 2.2). In 
all panels the solid blue circles correspond to sources that are de- 
tected in IRAC [3.6] (> 2cr), open squares are used for marginally 
detected sources (1-2<t), and gray open symbols are for IRAC- 
undctcctcd (< la) sources. The histograms in the right panels 
show the color distributions. Top: the 1775 — [3.6] color as a func- 
tion of magnitude in the [3.6] channel. The dashed line in the top 
right corner corresponds to the limiting magnitudes probed by our 
HST observations. The three symbols in the lower left show the 
typical uncertainties. The solid horizontal lines show the minimum 
and maximum color for a dust-free stellar population with a CSF 
history and ages between 10 Myr and the age of the universe at 
2 = 3.5. The dashed thick line is a fit to the detected sources and 
shows a trend of redder colors for brighter galaxies. The best fit 
corresponds to (4775 — [3.6]) = —0.25 X [3.6] + const. Middle: the 
UV-slope (f\ oc A' 3 ) determined from a linear fit to the available 
HST photometry of each sou rce (e.g., as done by IBouwens et all 
12012 : [Castcllano ct al. 2012]). This color also shows a trend of 
redder color s for brighter sources that is consistent with previous 
reports (e.g. Bouwcns ct al. 2012). The best fit line corresponds to 
= —0.19 X [3.6] + 2.87. In our stellar population modeling we 
assume that this t rend is caused exclu sively by a varying amount 
of dust reddening ( Mcurcr ct al. 1999 relation, Equation [TJ . Bot- 
tom : the 4775 — [3.6] color after dust reddening corrections are 
applied. These corrections are derived from the UV slope. The 
trend of redder colors for brighter sources is weaker, with a best fit 
corresponding to (1775 — [3. 6]) dust— corrected 

= -0. 10 x [3.6] + const. 
After this corrections, most sources lie within the range of colors 
expected for dust-free CSF models. The residual trend could be 
caused by a weak trend in age versus luminosity. 

maximum color (set by the age of the universe at z ~ 3.5) 
of a CSF model in the absence of dust reddening. The 
thick dashed line is a fit to the detected points only and 
shows a trend to redder UV-to-optical colors for brighter 
galaxies. 

The middle panel of the figure shows the UV-slope 
characterized by (3 (fx oc A* 3 ), as a function of the mag- 
nitude in the [3.6] channel. A similar trend to redder 
colors for brighter sources is observed. The UV-slope 
j3 was measured by fitting a power law to all the avail- 
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Fig. 2. — The observed colors of 2 ~ 5 and 6 LBGs. Left column 
(green points) show the colors of the V— dropouts (2 ~ 5). Right 
column (red symbols) shows the colors of the i— dropouts (2 ~ 6). 
Symbols are as in Figure [l] The same trends of bluer colors (both 
the UV-slope and the UV-to-optical colors) for fainter sources are 
present at these redshifts. 

able HST /ACS and WFC3/IR photometric points f or 
any given source (e.g., as done bv lBouwens et al.ll2012l or 
ICastellano et al.l|2012l). 

IBouwens et al.l (|2012f ) study the dependence of j3 
with UV luminosity and finds a similar tr end (see also 
iFinkelstein et all 120111: iMcLure et al.l 120111 ). As argued 
in that work. /3 can be affected by several factors but 
it is most strongly dependent on the total dust content 
(which sets the normalization applied to the dust curve 
assumed). 

The bottom panel of Figure [1] shows the result of ap- 
plying a dust reddening correction derived solely from 
the UV-slope to the UV-to-optical color. We find that 
the uncorrected trend in UV-to-optical colors versus [3.6] 
is best fit by the relation (i 775 - [3.6]) = -0.25 x [3.6] + 
const. After applying the dust correction the trend is 
now weaker (i 775 - [3.6])dust-corrected = -0.10 x [3.6] + 
const and most sources in the sample are in the color 
range that is covered by our dust-free CSF models. The 
trend in UV-to-optical colors (after removing the effect of 
dust) still shows a slight dependence on the [3.6] luminos- 
ity, but this is probably due to the fact that sources that 
are brighter in the rest-frame optical are on average older . 
This is similar to the results shown in lLabbe et all (|2007l) , 
who show that the UV and UV-to-optical color trends ex- 
hibited by star forming galaxies in the z ~ 0.7-3.5 range 
can be primarily explained by the effects of dust redden- 
ing, with a small contribution due to age. Figure [2] shows 
a very similar behavior at z ~ 5 and 6. 

In view of this result, and to alleviate the difficulties 
produced by the age-dust degeneracies, in the follow- 
ing modeling we establish the dust reddening directly 
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from t he observed UV slope (3 following the lMeurer et al.l 
<|1999D relation: 

A 1600 = 4.43 + 1.99/3 (1) 

where Aiqoo is the attenuation in magnitudes at 1600 A. 
The att enuations at o t her w avelengths are determined 
using a iCalzetti et al.l (|200Q[ ) extinction curve. Sources 
with j3 < —2.23 are assigned ^4i6oo = 0. 

3.2. The Constant Star Formation model 

We have used the FAST (jKriek et al.ll2009l ) SED fit- 
ting code to fit the observed rest-frame UV + optical 
SEDs of the galaxies in our sample wi th a set of synthetic 
stella r population models (we use the lBruzual &: Charlotl 
120031 models) . In this section we present the results ob- 
tained when a smooth constant star formation (CSF) 
rate history is assu med. We have assumed a mctallicity 
Z = 0.2 Z Q , and a iSalpeterl (| 19551) IMF with cutoffs at 
0.1 Mq and 100 Mq. The set of models only include the 
stellar continuum fluxes and does not consider the effects 
of emission lines. These are explored in later sections. 

Additionally, and as described in the previous section, 
we do not leave the reddening by dust as a free parameter 
in our modeling, but rather fix it based on the observed 
UV-slope (3 (Equation[T]). hi doing so, we greatly simplify 
the stellar population modeling of high-rcdshift galaxies, 
by eliminating the very significant degeneracy between 
dust and age. While the impact of this degeneracy on the 
derived stellar mass is more limited, this degeneracy can 
have a big impact on the SFRs. For example, at z ~ 4, a 
model with zero dust reddening can produce SFRs that 
are lower by 1.8 dex compared to a model with maximal 
dust reddening (minimum age). Meanwhile, the effect on 
the stellar mass is generally < 0.4 dex. 

The results of the SED fitting procedure are shown in 
Figure [3] The figure shows the SFR versus M sta r for the 
z ~ 4 (top, blue), z ~ 5 (middle, green), and z ~ 6 
(red) sources. The dashed lines show the minimum and 
maximum values of M sta r at a given SFR that our mod- 
els allow for. The distributions are very similar at all 
rcdshifts, suggesting little evolution of this relation, con- 
sistent with previous results that use similar a ssumptions 
(|Stark et alll200l IGonzalez et al.H2010ll20Tl . It should 
be noted, however, that in previous works, the relation 
presented generally corresponds to SFRs that have not 
been corrected for dust extinction. Figure [3] shows the in- 
trinsic SFRs, i.e., after correction for dust reddening and 
dimming. The stellar masses have been derived from the 
same models used to derive the SFR. The fact that the 
relations still show little evolution is a consequence of 
the UV and optical colors being similar, implying sim- 
ilar properties of the galaxies as a function of redshift , 
including dust correc tions (see also IGonzalez et alJl2012t 
iBouwens et"~aTll2012D . 

A small fraction of the sample shows ages that are ex- 
tremely young, crowding the 10 Myr dashed line. These 
sources are generally sources that are marginally de- 
tected or non-detected in the IRAC channels. The very 
blue UV-slopes of these sources dominate the fitting, 
driving them to young ages. It is possible that they are 
indeed very young but this is hard to establish with the 
depths of the IRAC data used here. The analysis in 
this paper focuses on the more massive sources that are 




Fig. 3. — The SFR versus M s t a r relation for for our samples at 
z ~ 4 (top, blue), z ~ 5 (middle, green), and z ~ 6 (botom, 
red). Only sources with reliable cleaned IRAC photometry are 
considered (see text Section 2.2). Stellar masses and SFRs were 
derived assuming a CSF history. The dust reddening in the model 
has been derived directly from the UV slopes f) of each sources 
according to Equation[T] In all panels, solid circles indicate sources 
that are detected (> 2<r) in the Spitzer /IRAC [3.6] channel and 
hence their masses can be more accurately estimated; open squares 
indicate marginal detections (1 — 2<r); and open gray circles indicate 
IRAC undetected galaxies. The dashed lines indicate the minimum 
and maximum M a t ar that a galaxy can reach assuming that it 
has been forming stars at a constant rate, for the minimum age 
included in our models (10 Myr) and for the age of the universe at 
z ~ 3.5,4.5, and 5.5 (from top to bottom). 

detected in Spitzer /IRAC, so this small fraction of the 
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of galaxies rises with time (e.g, 


Papovich et alj 12011b 


IReddv et al.H2012HSmit et alJl2012 


), something that has 



Fig. 4. — Comparison of the SFRs derived using different model assumptions. The top row shows the comparison of the SFRs derived 
assuming rising star formation histories and constant star formation histories. The bottom row shows the effect that correcting for optical 
emission lines have on the SFRs. Overall, we find that the derived SFRs are not very sensitive to the model assumptions. 

population should not bias our main results. 

In the following sections we explore the impact of two 
important ingredients that have not been taken into ac- 
count in our previous modeling. We use the results from 
the CSF model as our base for comparison. Our first 
major consideration is the impact of assuming a rising 
SFH (which reproduces better the evolution of the UV 
LF at z > 3) on our SFR and stellar mass estimates. 
Our second major consideration is to explore the effects 
of emission lines on the derived physical parameters. 



3.3. The Rising Star Formation Model 

Under the assumption of a smooth star formation his- 
tory, the UV luminosity o f galaxies is related to their SFR 
(e.g.. the Madau relation; iMadau et al.lll998t see also the 
S FR functions der ived from the UV luminosity functions 
- ISmit et al.H2012l ). The exact form of the relation does 
not change strongly between different smooth SFHs if 
the ages considered are older than ~ 100 Myr. Based 
on this relation, the previously presented model with a 
CSF history, makes a clear prediction for the evolution 
of the UV LF: it predicts a LF that does not change 
with redshift. This has been thoro ughly ruled out usin 
large samples of LBGs at z > 4 (e.g lBouwens et afll200 
l201lMlOesc"h et al]l2010h . 

The observed UV LF evolves with cosmic time show- 
ing a brighter characteristic magnitude at lower redshifts 
at least until z ~ 3. A model that better matches 
the observed evolution is one in which the average SFR 



also been predicted in several numerical simulations (e.g 
iFinlator et~alll20111 : Uaacks et alJl2012af >. 

By following sources at a constant cumulative num- 
ber densit y n(< Mjjy) = 2 x 10 ~ 4 Mpc~ 3 as a function 
of redshift. iPapovich et al.l (|2011[ ) derive a best- fit expo- 
nentially rising SFH of the form SFR cx e (t - to)/r , with 
r = 420 Myr (they find a slightly better fit using a linear 
model but the differences are not large). A similar anal- 
ysis at multiple values for the number density yields a 
similar result, with a best fit r ~ 500 Myr at all densities 
(Oesch et al. in preparation). The goal of the following 
analysis is to study the effects of such a model on the 
determination of the stellar masses and star formation 
rates. The details of the star formation history used are 
not important as long as it agrees (or at least agrees bet- 
ter than the CSF model) with the evolution of the UV 
LF. 

We have re-determined the SFRs and stellar masses 
for the galaxies in our sample using a RSF model as 
described before (SFR cx c( t_t °'/ T ) , with r = 500 Myr. 
The parameter to is free in our model. For the RSF and 
the CSF models, the SFR can be derived directly from 
the UV luminosities (if the dust extinction is known - 
here we assume it can be derived from the UV-slope). In 
fact, Figure [4] shows that, as expected, the SFRs derived 
from the CSF and the RSF model are on average the 
same with little scatter. Since the final goal is to study 



7 



z~4 



z~5 



z ~6 



11 F — I 1 1 1 — 

x 1 • IRAC detected ( >2tr) «^ 



10 - 
9 - 



□ IRAC marginal (lu-2tr) q 
O IRAC undetected 




7 8 9 10 

^csf/^g 



7 8 9 10 

^csf/^g 



7 8 9 10 

M CSF /M Q 



Fig. 5. — Comparisons of the stellar masses we derive using different model assumptions. The horizontal axis shows the stellar masses 
determined from our CSF model. The vertical axis for the upper set of panels shows the stellar masses derived from a model with 
exponentially rising star formation rate (SFR oc e( t-t o)/ T , where t — to is the age, and r = 500 Myr). This model is explicitly constructed 
to reproduce the observed evolution of the UV LF. Both Mqsf and Mrsf are derived with models with stellar continuum only. In all 
panels, the solid symbols represent sources that are detected (> 2a) in IRAC [3.6], open squares show marginal detections (l-2<r), and 
open gray circles show IRAC undetected sources. The left column shows the z ~ 4 sample (blue), the middle is the z ~ 5 (green), and 
the right column shows the z ~ 6 sample. When only the detected sources are considered, the mean AfcSF and the Mj^gp are consistent 
with each other, but with a scatter of ~ 0.25 dcx. The IRAC undetected and marginal detections show a slightly larger mass for the RSF 
determinations. The vertical axis for the lower set of panels shows the determinations of the stellar mass using a model that assumes strong 
optical emission lines which evolve with redshift (see Section 3.4). A CSF is assumed. The strength of the emission lines assumed can be 
seen in Figure [6] The effects on the stellar mass are larger at z ~ 6 where they cause an average decrease in the stellar mass of ~ 0.26 dcx 
for IRAC detected sources. This is expected as both IRAC channels present large contributions of flux from emission lines at this redshift. 



the sSFR in our sample, we now focus on the effects of 
the RSF on the stellar masses. 

As can be seen in the Figure [5] (top row) there is a sig- 
nificant scatter (~ 0.25 dex) in the comparison between 
the CSF and the RSF stellar masses. The mean value, 
however, docs not change significantly when the IRAC 
detected sources are considered (solid circles) . The IRAC 
undetected sources, meanwhile, show a slight bias to- 
wards larger RSF masses, but this could be attributed 
to the poor constraints in the rest-frame optical, which 
is very important to derive st ellar masses. This i s con- 
sistent with a recent report by iReddv et al.l (|2012[ ) , who 
also find that the average M sta r does not change between 
models with CSF and RSF for a sample of galaxies at 
z ~ 2-3. It is howeve r, inconsistent with the reports 
of IPapovich et all (|2011l ). who finds larger stellar masses 
when a RSF is assumed. In their analysis, however, the 
age is fixed with a formation redshift z = 11, artificially 
fixing the ages of the sources. 

There is one important observation regarding the RSF 
models. A significant fraction of the sources that we 
have modeled result with ages that correspond to the 



maximal allowed in our set of models. This is a result 
of the fairly red UV-to-optical colors exhibited by these 
galaxies (after dust reddening corrections). The mod- 
els considered here only include stellar continuum light, 
and as a consequence, these red colors can only be inter- 
preted as older ages. Two possible ways to alleviate this 
tension are: the presence of emission lines that would 
contribute to the flux mea sured with Spitzer /IRAC 
(|Schaerer fc de Barrosl [20ToT) ; or a non-smooth compo- 
nent to the SFH such that the SFH is rising on average 
but there is a characteristic duty cycle that make the col- 
ors redder at any given age (iLabbe et al.ll2010at see also 
recent simulation results by. e.g.. Uaacks 1TaTll2012bl) . 
Since many of the sources that present maximal age at 
z ~ 4 arc expected to be free of emission line contam- 
ination given their redshifts, the actual solution to the 
problem is probably a combination of both scenarios. 

In summary, both the SFRs and the stellar masses ob- 
tained from the RSF are on average the same as for the 
CSF model. We would therefore not expect the average 
sSFR to be substantially affected by whether one adopts 
a RSF or CSF history for the stellar population model- 
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Fig. 6. — The possible effect of strong rest-frame optical emission lines in the Spitzer /IRAC photometry. In both panels the blue curves 
show the corrections, in magnitudes, that need to be applied to the [3.6] photometry to remove the flux from emission lines. The green 
curve is the equivalent for the [4.5] photometry. In the left panel we conservati vely assume a constant Ha rest - frame EW of 100 A. To 
estimate the strength of the other lines we assume the flux ratios presented by Anders & Fritzc-v. Alvenslcbcn (2003), for a mctallicity 
Z = 0.2 Zq. We also assume a flat underlying continuum (in F v ) . In the right panel the Ha EW is estimated from the trends observed as a 
function of redshift for sources with logi (M s tar /Mq) = 10-10.5 IjFumagalli et aLll2012n , This extrapolation is based on recent HST grism 
observations and predicts Ha EWo ~ 300 A at z ~ 6. The light blue, light green, and light red, solid histograms (arbitrary normalization), 
show the redshift distribution of our sample. These corrections result in lower stellar masses but, in general, they do not affect the SFRs. 
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3.4. The possible impact of optical emission lines. 

We turn now to the effect of the possible contribution 
of emission lines to the rest-frame optical photometry. 
These lines are not expected to affect significantly the 
SFRs derived from SED fitting but they can, in principle, 
have a strong effect on the values derived for the stellar 
masses. 

It has been suggested that a (possibly large) fraction 
of the flux detected with Spitzer/IRAC from galaxies 
at z > 4.5 could be coming from nebular regions asso- 
ciated wjtli_sta2^Jcj^^^ of emission lines 
(e.g.. iSchaerer fc de Barrosll20l61 ). Unfortunately, it is 
not possible to directly observe these emission lines in 
z > 4 galaxies with spectroscopy using current facilities. 

Nonetheless, we can infer the approximate strength of 
these emission lin e s indi rectly. For example, in a recent 
study, iShim et ail (|201lD used a spectroscopic sample at 
3.8 < z < 5.0 to study the Ha EW from broad band 
photometry. In this redshift range, Ha falls in the IRAC 
[3.6] filter, adding to the rest-frame optical continuum 
flux at this wavelength. However, Ha does not contribute 
to light in the adjacent K-band or IRAC [4.5] filters, 
making it possible to infer the EW of this emission line. 
IShim et all (pOllh find a median Ha EW of - 480 A. 
This technique, however, is only sensitive to objects with 
large Ha EW (> 350 A). While this experiment shows 
that there is a fraction of the star-forming population 
at z ~ 4 with strong line emission, the objects selected 
by IShim et al.l (|2011[ ) are not expected to be typical for 
the full population. Independen t evidence f or ext reme 
line emitters at z > 2 come from lAtek et all (|201lD who 
identify a large sample of such objects at z ~ 2.5 based 



on HST WFC3/IR Grism data. 

Since the direct observation of emission lines at these 
rcdshifts is not currently possible, we need to make as- 
sumptions based on our understanding of the se lines at 
lower redshift. At z ~ 2 — 2.5, for example, lErb et all 
(|2006f ) find that the rest-frame EW of Ha for galaxies 
with log 10 (M star ) - 10.0-10.5 is - 100 A. The EW shows 
an increase with decreasing M s t ar . Moreover, this EW 
also seems to increase toward higher redshift. In a study 
base d on HST Grism spect r oscopy from the 3DHST sur- 
vey (jBrammer et al.ll2012f) . iFumagalli et al.l (|2012D find 
that the Ha EWs of 10 10 -10 10 5 M Q galaxies at z ~ 0-2.5 
are best fit by the following relation: 



EW{z) ~ 15.8 x [l + z) l - 52 k, 



(2) 



The typical strength of the optical emission lines at 
higher redshifts and lower masses will likely remain un- 
certain for some time. In the following, we estimate the 
effects that seem plausible based on these observations. 

We take two different approaches. In our first, more 
conservative approach, we will assume that the mean Ha 
EW observed at z ~ 2-2.5 (jErb et al.ll2006T) remains con- 
stant at 100 A at all redshifts z > 3 and stellar masses. 
Given the trends in the EW described earlier with stel- 
lar mass and redshift, this is likely an underestimate of 
the EW exhibited by real galaxies at very high-redshift. 
A more extreme approach is to assume that the trends 
observed at z ~ — 2.5 continue to higher rcdshifts fol- 
lowing the same extrapolation ( Equation^ . We call this 
the maximal emission line model. For each of the mod- 
els [conservative and maximal) we derive the strength 
of al l the other emission lines based on th e flux ratios 
from I Anders fc Fritze-v. Alv enslcbcnj (|2003f ) assuming a 
Z = 0.2 Zq mctallicity and a flat F v underlying opti- 
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cal continuum. Then we calculate the contribution of 
the emission lines to the broadband magnitudes assum- 
ing the redshifts derived from the CSF models in Section 
3.2. Figure [6] shows the emission line contributions (in 
magnitudes) that we estimate as a function of redshift. 
The solid histograms (arbitrary normalization) show the 
redshift distribution of our samples. 

In a previous study ((Gonzalez et al.l 120111) . we exam- 
ined the median SEDs of this sample by stacking the 
photometry in bins of UV luminosity. The SEDs of z > 4 
galaxies show a consistent excess in their [3.6] flux over 
their [4.5] flux which, as argued in that work, can be ex- 
plained by the effect of emission lines. It was shown in 
that work that a simple model with Ha EW = constant 
= 300 A, can simultaneously reproduce the observed UV- 
to-optical and [3.6] — [4.5] colors exhibited for the sam- 
ple. This EW is very similar to the value predicted by 
the maximal model at z ~ 6. 

Next, we subtract the emission line contribution from 
the photometry, and we refit the SEDs with CSF and 
RSF models. As expected, the SFRs derived in this 
way are unchanged in both cases (Figure^! bottom row). 
This is because the SFRs of star forming galaxies (with 
ages > 100 Myr) depend almost exclusively on the UV 
fluxes, which are unaffected by these emission lines. The 
effects on the stellar masses can be seen in Figure [5] (bot- 
tom row) . This figure only shows the effects for the max- 
imal model. The effects from the conservative model are 
much smaller. 

The impact of the emission lines on the stellar mass 
of galaxies is largest for the z ~ 6 sample. This is ex- 
pected, since at this redshift there are strong emission 
lines affecting both the [3.6] (H/3 and O III) and the [4.5] 
(affected by Ha) channels. The average change in stellar 
mass estimates AM 8 t a r for IRAC detected sources at this 
redshift is ~ 0.26 dex for the maximal model (and only 
0.1 dex for the conservative model). Correcting for emis- 
sion lines results in lo wer stellar masses and hig her sSFRs 
for galaxies (see also ICurtis-Lake et aLl l2~012). The im- 
pact of this on the evolution of the sSFR, however, is not 
straightforward to assess, since correcting for the effect 
of the emission lines on the mass shifts sources to lower 
mass where the sSFR also may be higher. 

4. THE SSFR AT Z > 2 

In the preceding sections we have studied the effects 
that different stellar population modeling assumptions 
have on the SFR and M s t ar estimates for high-rcdshift 
LBGs. We have shown that the SFRs do not depend 
on whether a rising SFH or CSF is considered in the 
modeling. The SFRs are also not significantly affected 
when emission lines, which affect the rest-frame opti- 
cal fluxes, are taken into consideration. For the simple 
smooth RSF that we consider, the average M sta r derived 
for these galaxies does not change systematically with 
respect to the CSF model (although there is a scatter of 
~ 0.25 dex in the values determined). When the emission 
lines are considered, however, the stellar masses could be 
~ 0.26 dex smaller at z ~ 6. 

Figure [7] shows the sSFR determined from our sam- 
ple as a function of redshift z > 4 (see also Table . 
The left panels shows the mean sSFRs derived from a 
CSF model and the right panels show the results using 
the exponentially rising SFH. The top panels correspond 



to stellar continuum only models and the bottom pan- 
els consider the effects of emission lines assuming the 
average contribution expected from the maximal model. 
The mean sSFRs were determined for galaxies in two dif- 
ferent M s tar bins. The blue circles and green hexagons 
correspond to galaxies with M sta r,csF ~ 5 x 10 9 M Q (the 
width of the bins are Alog(M star ) = ±0.3 dex). The or- 
ange pentagons and purple triangles show the sSFR for 
galaxies that are ~ x5 less massive. 

Galaxies with best fit ages < 10 Myr (the minimum 
in our grid of models) were rejected before estimating 
the mean values. The corresponding sSFR for the galax- 
ies that were rejected is sSFR ~ 100 Gyr -1 . The effect 
of keeping such galaxies in the more massive sample is 
shown by the gray symbols in Figure [7] The differences 
are only significant at z ~ 6. This may indicate that the 
z ~ 6 sample is dominated by very young sources (with 
high sSFRs), however, it is unlikely that such a sudden 
change in the characteristics of the population may arise 
in the period from z ~ 6 to z ~ 5 (a time span of only 
~ 300 Myr). This is especially questionable when con- 
sidering the small size of the z ~ 6 sample (< 10 sources 
per bin) and the decrease in the signal-to-noise of the 
data for these more distant galaxies (compared to the 
z ~ 4 and 5 samples). This redshift bin is also the most 
sensitive to the line emission model assumed. 

Comparing the mean sSFRs estimated for both mass 
bins, the data may suggest that the less massive galax- 
ies present slightly larger sSFRs. However, the differ- 
ences are small (< 0.1 dex) and very uncertain. For 
both mass bins and regardless of the modeling assump- 
tions presented here, the sSFR shows very little vari- 
ation with redshift fr om z ~ 4 to z ~ 5, consistent 
with previous reports (|Stark et al.ll2009l ; iGonzalez et al.l 
l2010HMcLure et al.ll2011l ). However, there is large scat- 
ter among the different z ~ 6 determinations. These 
variations reflect the limitations of the sample but also 
the fact that the sSFR at z ~ 6 is highly sensitive to the 
modeling assumptions. 

It should be noted that the galaxies that make up the 
1 x 10 9 Mq and 5 x 10 9 M Q samples are not the same 
galaxies in the top and bottom panels. In the top panel 
the Mgtar from CSF models with no emission lines is 
considered whereas in the bottom, the masses used cor- 
respond to the ones derived when the maximal model 
of emission lines is assumed. This is important because 
if the same galaxies were considered, then their sSFRs 
should be larger when the emission lines are considered 
(because their masses are lower). Most importantly, this 
effect would be stronger at z ~ 6, suggesting a stronger 
evolution of the sSFR with redshift. This does not seem 
to be the case when galaxies within the same mass bin 
are compared. 

The mean sSFR values derived with all the model as- 
sumptions are slightly larger but close to ~ 2 Gyr -1 , 
similar to previous reports, and consistent wi t h the 
values reported at z ~ 2-3 (]Daddi et al.l 120071 120091 ; 
iReddv et all 120121 see Figure El). A weighted best fit 
to the evolution of the sSFR as a function of redshift is: 
log(sSFR(z)/Gyr -1 ) = -0.1(±0.1)+0.9(±0.1) log(l+z). 
This fit w as obtained c onsid ering the lo w redshift re- 
sults from iDaddi et all lj2007l z ~ 2) and IReddv et all 
(|2012t z ~ 2-3) as well as the values derived here for 
the mean sSFR at z > 4 (taking the RSF with max- 
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Fig. 7. — The mean specific SFR as a function of redshift for our samples at z ~ 4, 5, and 6. The mean sSFR was determined 
for galaxies in two different mass bins 1 x 10 9 Mq and 5 X 10 9 M©. Top-left: the mean sSFR derived assuming a CSF history and 
models with stellar continuum only (i.e., ignoring the effects of emission lines). The mean values were estimated after removing extreme 
values (sSFR ~ 100 Gyr -1 ) which correspond to the minimum ages in our grid of models (10 Myr). This rejection only causes significant 
differences at z ~ 6, as can be seen by comparing the blue (with rejection) and gray circles (without rejection). Top-right: mean sSFR 
derived assuming an exponentially rising SFH and models with stellar continuum only. The mean sSFR values do not change significantly 
for the RSF model relative to that for a CSF model (see Top-left panel). Bottom- left: assuming our maximal model for the average 
strength of the emission lines, we re-model the galaxies first assuming a CSF SFH. The behavior is very similar to that shown in the top-left 
panel. The only differences are at z ~ 6, though they are not really significant within the errors. Bottom-right: the mean sSFR derived 
assuming an exponentially rising SFH and also assuming our maximal model for the average strength of the emission lines. Again very 
similar behavior is seen. Overall, as can be seen in the four panels, the mean sSFR at z ~ 4 and 5 varies weakly with both redshift and 
mass. The data may suggest that galaxies in the lower mass bin may have slightly larger mean sSFR values but the evidence is weak and 
the difference is at most 0.1 dex. The only significant differences are observed at z ~ 6. It should be noted, however, that the sSFRs at 
z 6 are the most poorly constrained given the small number of sources (generally about 8 per bin, see Table The z ~ 6 values are 
also the most sensitive to the different assumptions regarding the average strength of the emission lines, and whether extreme values are 
rejected from the mean or not. 



imal emission line model as our fiducial values). The 
weights used correspond to the random errors. Since 
the random errors are not available for many determi- 
nations in the literature, the weights were normalized 
based on the sizes of the samples. This effectively re- 
duces the weight of the z ~ 6 point significantly. Our 
result suggests a mean sSFR that increases slightly with 
redshift, however, the evolution observed is much weaker 
than that seen at z < 2. It is also inconsistent with 
the slope expected from theoretical models which gen- 
erally follow closely the specific halo a ccretion rate, re- 
sulting in sSFR(z) oc (1 + z) 2A (e.g ., IWeinmann et al.l 
[20111: iBouche et all 120101: lDavll2008h . If the CSF model 
without emission lines is assumed instead, the derived 
evolution is even weaker: sSFR(z) oc (1 + 2:)°- 6±01 

5. DISCUSSION 

Many of the earliest studies of the sSFR at high- 
redshift had suggest ed that the sSFR did not evolve 
strongly at z > 2 dGonzalez et al.l 120101 : iStark et al.l 
I200a iMcLure et all [20111 ) . This is very different from 



the behavior observed at lower redshifts and also from 
the expectations of theoretical models that consistently 
predict a sSFR that declines monotonically with cosmic 
time when halos of a constant mass are s tudied (e.g ., 
IWeinmann eF^I|[20Tl1: IBouche et alll2010l: [Da^[200l . 
Most earlier studies, however, did not determine the 
SFRs and M sta r m a way that was clearly self-consistent, 
did not adopt SFHs which are consistent with the evo- 
lution of the UV LF, and also did not account for the 
effects of nebular emission lines in the stellar population 
modeling. 

More recently, however, there have bee n some efforts to 
redres s these shortcomings. For example. iBouwens et al.l 
(|2012D focus on correcting previous sSFR determinations 
to reflect the latest estima tes of dust extinc t ion. A nother 
example is the work of Ide Barros et a D ([2TJT21 ) where 
there is an effort to correct for the impact of the emis- 
sion lines on the inferred stellar masses. In general, all 
these recent studies, including the present work, have 
suggested somewhat higher sSFRs at z > 4. 

It is worth noting that despite the growing consensus 
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TABLE 2 

Current Estimates of the mean Specific Star Formation Rate for 



4-6 Galaxies 







2 ~ 4 


2 ~ 5 




2 ~ 6 


Model 


iVbin 


login (sSFR/Gyr- 1 ) 


Nbin loe-mfsSFR/Gyr" 1 ) 


N hin 


loKmfsSFR/Gyr" 1 ) 








5 X 10 a Mq 






CSF no emission lines 


53 


0.49(±0.05) 


17 0.54(±0.09) 


7 


0.20(±0.06) 


CSF with emission lines 


49 


0.54(±0.05) 


16 0.62(±0.09) 


1 


0.35(±0.20) 


RSF no emission lines 


66 


0.48(±0.03) 


20 0.48(±0.04) 


9 


0.94(±0.22) 


RSF with emission lines 


63 


0.51(±0.03) 


20 0.52(±0.05) 


5 


1.10(±0.30) 






M s tar = 


1 x 10 9 M Q 






CSF no emission lines 


98 


0.59(±0.04) 


28 0.60(±0.07) 


5 


0.62(±0.18) 


CSF with emission lines 


98 


0.62(±0.05) 


31 0.66(±0.07) 


3 


0.64(±0.24) 


RSF no emission lines 


121 


0.65(±0.04) 


34 0.70(±0.08) 


11 


1.01(±0.18) 


RSF with emission lines 


122 


0.67(±0.04) 


37 0.79(±0.08) 


11 


1.09(±0.19) 



Note. — Mean values of the estimated log 10 (sSFR/Gyr _1 ) for our samples using different model assumptions. 
The sSFR was estimated for two stellar mass bins centered at log ln (M star /M0) = 9.7 and 9.0. The width of each 
bin is ±0.3 dex. The N\,\ n column indicates how many sources fall in each bin at each redshift for a given model. 
Extreme values of the sSFR (~ 100 Gyr -1 , which corresponding to minimum age models) were rejected before taking 
the mean. This latter choice only makes a difference at z ~ 6 (see Figure (TJ. 



that the SFRs at z > 4 increase, there is a difference in 
the magnit ude of evolution derive d. In particular, while 
studies like Ide Barros et al.l (|2012D find an order of mag- 
nitude larger sSFRs at z > 4 than is found at z ~ 2, the 
sSFRs we find at z > 4 are only larger than the z ~ 2 
values by factors of ~ 1.5-2.0. What are the reasons for 
these differences? Largely, this is the result of the differ- 
ent assumptions various authors use in modeling the ob- 
served photometry. Differences in assumptions can have 
an impact on the derived sSFRs due to degeneracies be- 
tween different model parameters. While we adopt the 
simplifying assumption that the observed UV continuum 
slope can be used to estimate the dust ex tinction in high- 
rcdshift galaxies, Ide Barros et al.l (|2012| ) do not impose 
any constraints on the dust, metallicity, mass, or age in 
modeling the photometry of high-rcdshift galaxies. By 
not restricting the model par ameter space th r ough vari- 
ous simplifying assumptions, Ide Barros et al.l (|2012f) ob- 
serve considerable scatter in the properties of many of 
their sources and also find a large fraction of sources 
to have very young ages (a result wh ich we consider 
proba bly unphysical and may result in Ide Barros et al.l 
120121 reporting very high sSFRs). On the other hand, 
with our approach, we use the UV continuum slope to 
set the dust extinction, effectively breaking many of the 
model degeneracies. While one can debate which ap- 
proach yields the most accurate results, we prefer our 
approach due to the strong evidence at both z ~ and 
z ~ 2 that the UV continuum slope c orrelates on ayerage 
with the observed dust extinction dMeurer etaLl [19991: 
Burgarella et~al] [20051: lOverzier et al.ll2011t iDaddi et al.l 
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■garena et ai-iizuuat iu 
tt lReddv et al.ll2012jh 



Despite the change in the average stellar masses when 
the emission lines are considered, the sSFR at a con- 
stant stellar mass does not change especially strongly 
with cosmic time. With or without the consideration of 
emission lines, the sSFR shows only a weak evolution at 
z > 2 (Figures [7] and [8} . A fit to the maximal emission 
line estimates of the sSFR, which inc ludes the values at 
z ~ 2 — 3 reporte d in the literature (jDaddi et al.l 120071 : 
iReddv et al.ll2012| ). yields a best-fit to the evolution of 
the sSFRoc (l + z)°- 9±01 , indicating that the sSFR is 
higher by ~ 2.1 x at z ~ 6 than at z ~ 2. The ex- 



ponent of this relation is lower if a correction is not 
made for the emission lines (0.6±0.1 for the CSF model). 
Still, the evolution at z > 2 is less than expected from 
the lower-z trends and inconsiste nt with the predictions 
from numerical simulat ions (e.g., iNeistein fc Deke3l2008l : 
IWeinmann et al.ll201l[ ). 

6. SUMMARY 

We take advantage of the ultra-deep and wide-area 
ACS + WFC3/IR + Spitzer observations of the HUDF 
and ERS fields to derive flux measurements for a size- 
able sample of z ~ 4-6 galaxies and use these measure- 
ments to more thoroughly quantify their SFRs, stellar 
masses, and other properties. We have explored the ef- 
fects of using a reddenin g law based on the UV colors 
only ()Meurer et al.l I1999T) and investigated the impact 
of the star formation history by alternatively consider- 
ing constant star formation (CSF) and smoothly rising 
star formation (RSF) models. We have also studied the 
effects of optical emission lines on the SFR and stellar 
masses assuming two different models for the emission 
line strengths. Our goal has been to quantify the impact 
of these assumptions in estimating the specific star for- 
mation rate of LBGs which has pre viously been reported 
to show little evo l ution at z ~ 2-7 (|Gonzalez et aLll2010l : 
IStark et al.l[2009t IMcLure et al.H 2011h. The RSF history 
explored here matches the evolution of the UV LF (e.g., 
iPapovich etaLl 120111: iSmit et al.l I2012D and is in good 
agreement with the predictions fro m recent smooth par - 
ticle hydrodynamics simulations (jFinlator et al.l 120 111 ) . 
Our main findings arc the following: 

• At a given redshift, LBGs that are fainter in the 
rest-frame optical show bluer UV colors. This is 
consistent with the previously rep orted trends of 
bluer UV co l ors for fainter galaxie s (jBouwens et al.l 
[2009L [20121 IWilkins et alllMl . The UV-to- 
optical colors are al s o bluer for fainter sources 
(|Papovich et al.l 120041: iGonzalez et al.ll2012t l. Such 
luminosity-dependent trends could arise from a 
dependence of the dust content on mass (or lu- 
minosity) or from reddening caused by aging 
with a smooth SFH (see also lLabbe et al.l 120071: 
iBouwens et al.l [200l 120121: IGonzalez et al.ll2012D . 
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Fig. 8. — The mean specific SFR as a function of redshift for galaxies with estimated stellar masses log-^Mstar/A^Q) = 9.4—10, corre- 
sponding to our 5 X 10 9 Mq bin. The values at z < 4 are taken from the literature (Damen et al. 2009; Noeske et al. 2007; Daddi et al. 2007 



spoi 

I Reddy et al .11201 21) . The solid gray points correspond to previous reports from the literature at z > 4 l lStark"* eta! . 2009; Gonzale z et al.l 
120101 ) and the open gray points correspond to these values after a simple correction for dust extinction is made (Bouwens et al. 2012]). The 
values marked by the open hexagons at z > 4 are our new results from our stellar population modeling of the full rest-frame UV and optical 
SED. These results were derived assuming an exponentially rising SFH (RSF) and the average emiss i on lin e flux is subtracted based on the 
maximal model described in Section 3.4. The dust reddening was derived using the Mcurer et al. (1999) relation and the measured UV 
slopes, /9. The error bars correspond to random errors. Our estimate of the systematic error, is given as well in the upper left. The mean 
sSFR depends weakly on the modeling assumptions. The largest fluctuations are observed at z ~ 6 but this data point is very uncertain 
given the small sample, the lower quality of the data, and the large impact that emission lines can have at this redshift. A weighted best fit 
to the ev olution of the mea n sSFR at z > 2 is: log(sSFR( z)) = — 0.1(±0.1) + 0.9(±0.1) log(l + z) (which incorporates the values derived at 
z ~ 2 bv lDaddi et~aTH2(307l . z ~ 2-3 bv lReddv etS]|20ia and the z > 4 from this work). The derived dependence of the sSFR on redshift 
is much weaker than that observed at z < 2 and that expected from the theoretical expectations at z > 2. 



Assuming that the UV slope reddening is caused 
by dust alone a s in Equation [T1( Mc urer et alJll99& 
and following a lCalzetti et al.lbOQOt dust law), the 
dust-corrected UV and UV-to-optical colors can be 
reproduced by dust-free CSF models. A residual 
trend towards bluer colors at lower luminosities can 
be cause d by the effects of agi ng of the population 
(see also iGonzalez et aT1l2012[ ). 

We find that the SFRs derived assuming a CSF 
history do not change when a smoothly rising SFH 
(RSF) is used instead. The RSF that we use 
is chosen to match the evolution of the UV LF 
(SFR oc e(*- to) / T ; r = 500 Myr). The mean stel- 
lar masses are also unchanged when the RSF is 
assumed, although there is a significant scatter of 
~ 0.25 dex in the M starj csF - M star ,RSF relation. 

We explore the effects of emission lines assuming 
two different models for the strength of the lines 
as a function of redshift. Regarding the SFRs and 
stellar masses, the impact of adding emission lines 
with strengths consistent with the trends at lower 
rcdshifts is important, but modest. The inclusion 
of emission lines only affects our derived stellar 
masses. For the maximal emission lines model, the 
stellar masses are up to 0. 26 dex lower at z ~ 6 (an d 
the sSFR larger; see also iCurtis-Lake et al.ll2012l ). 
For the model with lower emission line strengths, 
the masses change by only ~ 0.1 dex. 



• We estimate the sSFR for galaxies in two stel- 
lar mass bins in our samples at z ~ 4, 5, and 
6. There is a slight hint that the sSFRs of lower 
mass sources may be slightly larger than those of 
more massive systems. In particular, for the CSF 
model, 1 x 10 9 M & galaxies have sSFRs that are 
~ 0.1 dex larger than sources with 5 x 10 9 M (cf. 
IBouwens et al.l [2012D . However, the uncertainties 
are still large, and therefore the current results are 
also quite consistent with no change in the sSFR 
with mass. Using a RSF model does not change 
this results. 

• At a fixed stellar mass, the sSFRs we derive only 
show a modest amount of evolution over the red- 
shift range z ~ 6 to z ~ 2. Emission lines lower 
the derived stellar masses and increase the sSFRs 
of individual galaxies, but for an ensemble of galax- 
ies at constant stellar mass, the sSFR only shows 
a modest degree of evolution from z ~ 6 to z ~ 2. 
This is somewhat in co ntrast to the conclusions by 
Ide Barros et al.l (|2012h who argue for stronger evo- 
lution. The sSFR at z > 2 shows a dependence on 
redshift sSFR oc (1 + z)°- 9±01 , i.e., a factor ~ 2.1 x 
higher at z ~ 6 than at z ~ 2 (assuming the 
RSF model with maximal emission line strength). 
This evolution is much slower than observed at 
z < 2 and would appear to be inconsistent with 
the expectations from simulations, which predict a 
steeper monotonically-decrcasing slope of the sSFR 
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with redshift at z > 2 (sSFR oc (1 + z) 2 - 4 ) 
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